Magnetic properties of the 5 = 1/2 distorted 
diamond chain at T = 

Kiyomi Okamotof , Takashi Tonegawa$ and Makoto Kaburagi§ 

f Department of Physics, Tokyo Institute of Technology, Oh-Okayama, Meguro-ku, 
Tokyo 152-8551, Japan 

I Department of Mechanical Engineering, Fukui University of Technology, Gakuen 
3-Chome, Fukui-Shi 910-8505, Japan 

§Faculty of Cross-Cultural Studies, Kobe University, Tsurukabuto, Nada-ku, Kobe 
657-8501, Japan 

Abstract. We explore, at T = 0, the magnetic properties of the S = 1/2 
antiferromagnetic distorted diamond chain described by the Hamiltonian 7i = 

J2jLl {Jl ' S>3j + Ssj ■ Szj+i) + J^Ssj+l • &3j+2 + ^3 (&3j-2 1 &3j + Saj 1 ^+2)} 

-HY4L1 S i with Jl ' J2 ' J 3 - °> wnich wel1 models A 3 Cu 3 (P04)4 with A = Ca, Sr, 
Bi 4 Cu3V 2 0i4 and azurite Cu 3 (OH) 2 (C03)2- We employ the physical consideration, 
the degenerate perturbation theory, the level spectroscopy analysis of the numerical 
diagonalization data obtained by the Lanczos method and also the density matrix 
renormalization group (DMRG) method. We investigate the mechanisms of the mag- 
netization plateaux at M = M s /3 and M = (2/3)M s , and also show the precise phase 
diagrams on the (J2/J1, J3/J1) plane concerning with these magnetization plateaux, 
where M = &i an d Ms is the saturation magnetization. We also calculate the 

magnetization curves and the magnetization phase diagrams by means of the DMRG 
method. 



PACS numbers: 75.10.Jm, 75.40. Cx, 75.50.Ee, 75.50.Gg 
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1. Introduction 



Low-dimensional quantum spin systems have attracted increasing attention in recent 
years. A few years ago, Ishii, Tanaka, Mori, Uekusa, Ohashi, Tatani, Narumi and 
Kindo pQ have reported the experimental results on a trimerized S = 1/2 quantum spin 
chain Cu3Cl6(H 2 0)2 • 2H8C4SO2. From the structure-analysis experiment they have 
proposed a model shown in figure 1 for this substance pQ. The Hamiltonian of this 
model is written as 



Ti = Tio + Jii 
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where Si is the S = 1/2 operator at the Zth site and N is the total number of spins 
in the system. All the coupling constants, J\, J 2 and J3, are supposed to be positive 
(antiferromagnetic). Hereafter we set J 2 = J2I J\ an d J 3 = J3/ J\- The magnetic field is 
denoted by H. We can assume J\ > J3 without loss of generality because the J\ < J3 
case is equivalent to the J\ > J3 case by interchanging the role of J\ and J3. 




Figure 1. Sketch of the distorted diamond (DD) chain. Solid lines denote the coupling 
Ji, wavy lines the coupling J2, and dotted lines the coupling J3. We can assume J\ > J3 
without loss of generality. 



This model has been named the "distorted diamond (DD) chain model" [2]. The 
name "diamond" comes from the mark of the playing cards, and "distorted" from the 
fact Ji 7^ J3 in general. The symmetric case (Ji = J3) has first been proposed by 
Takano, Kubo and Sakamoto (TKS) .3] earlier than Ishii et al. jTj and ourselves 
We have investigated the ground state properties of the 5=1/2 DD chain model 
when H = by an analytical method as well as the level spectroscopy (LS) analysis 
of the numerical diagonalization data obtained by the Lanczos technique. Later the 
theoretical studies on the 5=1/2 DD chain have also been reported by ourselves [HE], 
Sano and Takano [0] and Honecker and Lauchli 
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Unfortunately the substance which was thought to be Cu 3 Cl 6 (H20) 2 • 2H 8 C 4 S02 at 
first P has been proved to be Cu 2 CLrH 8 C4S0 2 , the lattice of which is the two-leg zig-zag 
chain with the bond-alternation by the same group Now, however, the 5=1/2 DD 
chain model is again in spotlight as a model for A 3 Cu 3 (P04)4 with A = Ca, Sr Plll)|lll| . 
Bi 4 Cu3V 2 0i4 [121 and azurite Cu 3 (OH) 2 (C0 3 ) 2 . 

For A 3 Cu 3 (P0 4 ) 4 with A = Ca, Sr, Drillon et al. [3110] measured the magnetic 
susceptibility x(T), the specific heat C{T) and the magnetization process M(H), and 
proposed the model without J 2 interactions. Ajiro et al. ^1] performed the high field 
magnetization process measurement up to 40T and the neutron diffraction experiment, 
and proposed the model with J 2 interactions. In these reports, a wide magnetization 
plateau at M — M s /3 was observed, where M is the z component of the total spin, 
defined by M = J2iL± Sf, and M s is the saturation magnetization. The behavior of x{T) 
and neutron diffraction results suggest that these substances are ferrimagnetic above 
the three-dimensional ordering temperature. 

Sakurai et al. [E] reported x(T), M(H) (up to 28T), C(T) and the 51 V NMR of 
Bi 4 Cu 3 V 2 0i 4 . Unfortunately, they could not obtain any conclusion on the existence of 
the M = M s /3 plateau because M ~ 0.27M S even when H = 28 T. The ferrimagnetic 
behavior was not seen in x(T) above the three-dimensional ordering temperature. 
It seems that this substance has the spin-fluid (SF) ground state above the three- 
dimensional ordering temperature. 

Very recently, Kikuchi and co-workers have reported the experimental 

results on the magnetic and thermal properties of azurite Cu 3 (OH) 2 (C0 3 ) 2 . Following 
their reports, azurite has the spin-fluid (SF) ground state and the wide magnetization 
plateau at M — M s /3. 

Thus, above the three-dimensional ordering temperature, the ground states of 
A 3 Cu 3 (P0 4 ) 4 with A = Ca, Sr seem to be ferrimagnetic, whereas those of Bi 4 Cu 3 V 2 0i 4 
and azurite Cu 3 (OH) 2 (C0 3 ) 2 seem to be SF. In view of these situations, we think it is 
very important to investigate the magnetic properties of the S = 1/2 DD chain model 
in more detail. Throughout this paper we consider the T = case. 

This paper is organized as follows. The ground-state properties at zero magnetic 
field are reviewed in §2. The magnetization plateaux at M — M s /3 and M = (2/3) M s 
are discussed in §3 and §4, respectively. The numerical results for the magnetization 
curves and the magnetization phase diagram are presented in §5. The last section §6 is 
devoted to discussion and summary. 

2. Review of the ground-state properties at zero magnetic field 

Discussions presented in this section are confined to the case where H = 0. We j2] have 
investigated the ground state of the S =1/2 DD chain model by an analytical method 
as well as the LS analysis of the numerical diagonalization data obtained by the Lanczos 
technique, and obtained the phase diagram shown in figure El There are three phases in 
this phase diagram; the ferrimagnetic (FRI) phase, the dimer (D) phase and the spin- 
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fluid (SF) phase. The magnitude S tot of the total spin S tot , defined by S tot = Y4L1 &i 
and S% ot = S to t{S to t + 1), is S to t = N/6 in the FRI phase, whereas S tot = in the D and 
SF phases. In the D phase, there is a finite energy gap between the doubly degenerate 
ground state and the first excited state, while in the SF phase there is no energy gap. 




Figure 2. Zero field ground state phase diagram of the 5 = 1/2 DD chain model. 
The J3 = 1 case is reduced to the model of TKS There are three phase in this 
phase diagram; the ferrimagnetic (FRI) phase, the dimer (D) phase and the spin-fluid 
(SF) phase. 

When J 2 = 0, we can readily know that the ground state is the FRI state with 
Stot = N/6 by use of the Lieb-Mattis theorem [TSJ. The physical pictures of the FRI 
state are shown in figure 01 where the picture (a) holds good for the J\ ~ J3 case, 
whereas (b) for the J\ ^> J 2, J3 case. The ellipses in (b) denote the effective 3-spin 
cluster which is formed when the interaction J\ is much larger than the interactions J2 
and J3. We represent the Sf = 1/2 and Sf = —1/2 states of the /th single spin by f« 
and li, respectively. Then, the ground state wave functions of the jth cluster in the 
case of J 2 = J 3 = are expressed as 

= 7/= (I T3j-lt3j!3j+l) - 2| T3j-l|3jt3j+l) + I l3j-lT3jT3j+l)) (4) 

and 

4>2,j = (I Uj-iUjhj+i) - 2 | Uj-ihjUj+i) + I hj-iUjUj+i)) (5) 

for Stot =1/2 and S^ 2 = —1/2, respectively, where S t ot is the z component of the 
total spin of the 3-spin cluster. When all the 3-spin clusters are effectively in the 4>i ( or 
02) state, the FRI state shown by figure 3(b) is realized. This physical consideration 
has been developed in jU[7j by use of the degenerate perturbation theory [J7J around 
the point J 2 = J3 = 0. In fact, the quantum phase transition between the FRI and SF 
phases takes place at J3 = J2 near (J 2 , J3) = (0,0), and it is of the first order. This 
shows very good agreement with the numerical results (see figure |5J). The quantum 
phase transition between the FRI and D phases is also of the first order. 
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(b) 



Figure 3. Physical pictures of the FPJ ground state for the S = 1/2 DD chain model. 
Open ellipses with open triangles in (b) denote the 3-spin clusters in the state fa. 



In the case of J 3 = 0, our DD chain is reduced to the J\-J\-Ji trimerized 
antiferromagnetic chain, the ground state of which is the SF state. In particular, in 
the case of Ji = J 2 , it is further simplified to the uniform antiferromagnetic chain. 





Figure 4. Physical pictures of the D ground state for the S — 1/2 DD chain model. 
Two spins in a shadowed ellipse effectively form a singlet dimer pair. 



The D phase is caused by the frustration, and is doubly degenerate due to the 
spontaneous breaking of the translational symmetry (SBTS), as shown in figure 0] The 
D phase is realized in the region where the frustration effect seems to be severe. We 
note that there is no frustration when J 2 = or J 3 = 0. The quantum phase transition 
between the SF and D phases is essentially the same as that in the 5 = 1/2 chain with 
next-nearest-neighbor interactions [T%l ll9|. and is of the Berezinskii-Kosterlitz-Thouless 
(BKT) type jZOlIZI], as has been fully discussed in our previous paper [2]. 

As we have mentioned in §1, the symmetric case J\ = J3, which is shown in figure 
Efa), has first been proposed by TKS [3]. They have found that the ground state 
is the FRI state, the tetramer-dimer (TD) state or the dimer-monomer (DM) state, 
depending on whether J 2 < 0.909, 0.909 < J 2 < 2 or 2 < J 2 . The value 0.909 of the 
FRI-TD boundary has been calculated by the numerical method, whereas the value 2 
of the TD-DM boundary is the exact one obtained by an analytical method. In the 
TD state, as shown in figure Efb), four spins of the diamond unit S 3 j- 3 , S 3 j-2, S^j-i 
and S 3 j (or S 3 j, S 3 j + i, S 3 j + 2 and S 3 j +3 ) form a tetramer, and two spins S 3 j + \ and 
&3j+2 (or S , 3 j +4 and S 3 j +5 ) form a singlet dimer. Obviously, the TD state is doubly 
degenerate due to the SBTS. In the DM state, on the other hand, two spins S^j+i and 
S 3 j + 2, connected by the J 2 interactions, form a singlet dimer and the remaining joint 
spins Ssj are completely free, as shown in figure Efc) . Thus, the DM state is 2 7V//3 -fold 
degenerate. It is very interesting that the ground state can be exactly expressed as the 
direct product of the local states in a finite range of the parameter region J 2 > 0.909. 

The TD state is very peculiar to the J\ = J3 case. Namely, if J\ > J3, the tetramer 
which consists of S 3 j_ 3 , -S , 3 j_ 2 , #3j-i and S 3 j is effectively decomposed into two dimers, 
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Figure 5. (a) The model of Takano, Kubo and Sakamoto (b) The tetramer-dimer 
(TD) state, where the rectangles represent tetramers and the ellipses singlet dimers. 
(c) The dimer-monomer (DM) state. 



the (Szj-3, Ssj-2) pair and the (S 3j -_i, Sy) pair. Two spins in each pair are connected 
by the J\ interaction which is larger than the J3 one. This decomposition leads to the 
doubly degenerate D state shown in figure |H Thus, the TD state is the special case of 
the doubly degenerate D state. 

The 2 Ar ' /3 -fold degeneracy of the DM state is lifted when J\ ^ J 3 . In other words, the 
effective interaction between the free spins S-^ and £3(7+1) appears through the dimer 
located between them. Honecker and Lauchli [7j have derived the effective interaction 
between S^j and Sky+i) as 



J, 



eff 



(Ji- h? n 



(6) 



J 2 V2 

which leads to the SF ground state. The effective interaction J eff vanishes when J\ = J 3 , 
which is consistent with TKS's result [H] on the DM state. Then, the DM state is also 
peculiar to the J\ = J 3 case and is the special case of the SF state. 



3. Magnetization plateau at M — M s /3 

In this section, we discuss the magnetization plateau at M — M s /3. Since the 
Hamiltonian has the trimer nature, as is shown in figure |B1 the mechanism for the 
M = M s /3 plateau is just the same as that for the trimerized chain (Figure E| 
investigated by Okamoto and Kitazawa [22J and by Honecker |23j. The M = M s /3 
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plateau can be realized without any SBTS. This is consistent with the necessary 
condition for the magnetization plateau by Oshikawa, Yamanaka and Affleck (OYA) [24J, 

n(S — (m)) = integer (7) 

where n is the periodicity of the wave function of the plateau state, S the magnitude of 
spins and (m) the average magnetization per one spin in the plateau. 



Figure 6. Single chain representation of the DD chain model. 



J J' J' J J' 

o o — o — o o — o 



Figure 7. Antiferromagnetic chain with the trimerization. 




Figure 8. Two mechanisms for the M = M s /3 plateau for the S = 1/2 trimerized 
chain: (a) the plateau A and (b) the plateau B. Open ellipses denote the trimers with 
S^t* = 1/2 and shaded ellipses denote the singlet dimers 

Here we explain the mechanism for the M = M s /3 plateau of the 5 = 1/2 trimerized 
chain. When < J <C J', each of three spins connected by J' forms an effective trimer 
with S£l z = 1/2 (or -1/2), the wave function of which is 0i (or 2 )- When the 
magnetic field is applied, all the trimers belong to the Sfot =1/2 state, which brings 
about M = M s / 3, as shown in figure a). The plateau due to this mechanism is named 
the "plateau A" . The mechanism of the plateau A state is essentially the same as that 
of the FRI state for the J\ 3> J2, J3 case when H = 0, which is already explained in §2. 
When J ^> J' > 0, on the other hand, each pair of two spins connected by J forms an 
effective singlet dimer pair and the remaining spins are nearly free. When the magnetic 
field is applied, the nearly free spins turn to the direction of the field, resulting also in 
M = M s /3, as shown in figure |Hfb). We call the plateau due to this mechanism the 
"plateau B". Thus, there are two mechanisms for the M = M s /3 plateau. It is easily 
seen that the change of the plateau mechanism occurs at J = J', where the trimerized 
chain is reduced to the uniform chain having no gap (i.e., no plateau) in the excitation 



Magnetic properties of the S = 1/2 distorted diamond chain at T = 8 

spectrum. This situation is most simply seen in the trimerized XY chain (2H| which is 
exactly solvable. 




Figure 9. Two mechanisms for the M = M s /3 plateau for the S — 1/2 DD chain: 
(a) the plateau A and (b) the plateau B. Shaded ellipses denote the singlet dimers and 



f 3) z 

open ellipses denote the trimers with 5 tot =1/2 



A similar situation is expected to hold for the present 5 = 1/2 DD chain model. 
The pictures of the plateaux A and B are depicted in figure My a) and (b), respectively, by 
use of the diamond form. When J 3 = 0, our model is reduced to the simple trimerized 
chain discussed above. Then, the change of the plateau mechanism occurs at J\ = J 2 
when J3 = 0. To find the phase boundary between the plateaux A and B for the general 
case, we have to rely on the numerical method. One of the most powerful methods is 
the level spectroscopy (LS) [2HII2ZI- Kitazawa [2H] has developed the LS for this kind of 
transition (belonging to the Gaussian universality class) by use of the twisted boundary 
condition (TBC) instead of the periodic boundary condition (PBC). This method has 
been applied to several magnetization plateau problems in a successful way |22j29p30|l3l] . 

In finite size systems with the PBC, the M = M s /3 plateau state changes smoothly 
from the plateau A state to the plateau B state when the parameter J2/J1 runs from 
to 00, although the plateau width takes the smallest value at the transition point. 
Under the TBC 

on the other hand, the plateau A state and the plateau B state can be easily distinguished 
even in finite size systems by the eigenvalue P of the space inversion operator 



(9) 



Sn-3+1- 

Since it is believed that the boundary condition does not affect any physical quantities 
in the N — > 00 limit, we distinguish the ground state by comparing the lowest energies 
with P = ±1. In other words, we can find the phase boundary between the plateaux A 
and B from the crossing of the lowest energies with P = ±1 as functions of the quantum 
parameters in the Hamiltonian. This is the physical interpretation of Kitazawa's TBC 
method |28j . He j2H] has consolidated the foundations of his TBC method by use of the 
bosonization, the renormalization group method and the conformal field theory. 

Figure ITU1 shows the crossing of the lowest energies with P = ±1 when N = 18 and 
J3 = 0.55. We can see that two curves cross with each other at j| (18) ~ 1.193, close 
to which the quantum transition point between the plateau A and B phases for the 
iV —y 00 system is expected to be located. An accurate value of ma y be estimated 
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Figure 10. Crossing of the lowest energies in the M = M s /3 sector with P = 1 (closed 
circles) and P = — 1 (open circles) for the N = 18 and J3 = 0.55 case under the TBC, 
when J2 is run. The crossing point of these two curves locates at J2 ~ 1.193. 



by extrapolating (N) to the iV 
assumed an iV dependence of J2 (N) as 



00 limit. Performing this extrapolation we have 



jt\N) = Jr + c 1 /N 2 + c 2 /N 4 



?(cr) 



(10) 



where C\ and c 2 are numerical constants, and have employed the finite size values, 
jf r) (12), J^ r) (18) and jf r) (24). Repeating this procedure with sweeping J 3 , we have 
obtained the M = M s /3 plateau phase diagram shown in figure ITT1 The estimated value 
of errors of the phase boundary is, for instance, J 2 = 1-193 ± 0.001 (or better) for 
J 3 = 0.55. 

We have performed the density matrix renormalization group (DMRG) calculation 
j221ES] to obtain the ground state magnetization curves and the magnetization phase 
diagrams for our S = 1/2 DD chain model, the details of which will be described in 
§5. Here, we present in figure H21 the width (see equation (29) below) of the 

M = M s /3 plateau for J 3 = 0.55 in the iV — » 00 limit; in this figure is plotted as 

a function of J 2 . The point where W(\/3) vanishes is the transition point (J 2 — 1.193), the 
result shown in figure IT2*1 being in very good agreement with the above-mentioned result 
obtained by the LS method with the TBC. We see a cusp in the curve at J 2 ~ 0.561. 
This cusp corresponds to the FRI-SF transition in figure El (see also figure 19(b) below). 
In fact, an infinitesimal small field brings about the M = M s /3 magnetization state in 
the FRI region, whereas a finite field is required for the M = M s /3 magnetization state 
in the SF and D regions. 

Let us now discuss the expectation value (S*) of each spin in the M = M s /3 plateau 
state. In the plateau A region, the approximate wave function of the jth 3-spin cluster 
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Figure 11. M = M s /3 plateau phase diagram of the 5 = 1/2 DD chain model. 




Figure 12. Normalized M — M s /3 plateau width Wt\is\/H a as a function of J2 when 
J3 = 0.55, where H s is the saturation field. The plateau width Wn/3) vanishes at 
J3 ~ 1.193, which well agrees with the transition point between the plateaux A and B 
in figure fTTI The cusp at J3 ~ 0.561 corresponds to the FRI-SF transition in figure|21 



IS 

= ~/= (I T3j-lT3iJ.3i+l) - 2| T3j-l|3jt3j+l) + I Uj-lhjUj+l)) (11) 

which leads to 

<S*-i> = (^+i> = \ (Ssj) = ~ (12) 
Thus, the expectation values (Sf) should be 

3' ~6' 3 J ' U' ~6' 3 J ' \3' ~6' 3/ ' 7 ( } 
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Figure 13. Behavior of (Sf ) of the N — 96 system (a) in the plateau A region, where 
{J2, J3) = (0.80, 0.55), and (b) in the plateau B region, where (J2, J3) = (1.50, 0.55). 

where three spins in a parenthesis form a 3-spin cluster. In the plateau B region, on the 
other hand, it is very easy to see that (Sf) should be 

{•••,(0,0), ^,(0,0), i (0,0), i (0,0), i ■••} (14) 

where two spins in a parenthesis form a singlet dimer. Figure IT31 shows the DMRG result 
for (Sf) for the N = 96 system in the plateau A region, where (J2, J3) = (0.80, 0.55), 
and that in the plateau B region, where (J2, J3) = (1.50, 0.55). Then, the behavior 
of (Sf) in this figure is quite consistent with our picture of the mechanisms for the 
M = M s /3 plateau. 

4. Magnetization plateau at M = (2/3)M s 

Our incipient discussions on the M = (2/3)M s plateau of the 5 = 1/2 DD chain model 
are given in a previous paper here we aim at discussing it in much more detail. 

Since the unit cell of the present model consists of three S — 1/2 spins, the OYA 
condition j21] tells us that if the SBTS does not occur, the possible magnetization 
plateau for < M < M s is only at M = M s /3. Thus, the SBTS is necessary for the 
realization of the M = (2/3) M g plateau. In the 3-spin cluster limit where J2 = J3 = 0, 
the M = (2/3) M s magnetization state is realized when a half of the 3-spin clusters 
are in the S^l z = 1/2 state and the remaining half are in the S^ z = 3/2 state. The 
configuration of these two states are completely free when J 2 = J 3 = 0, resulting in the 
2 Ar / 3 -fold de generacy. This high degeneracy is lifted when we introduce J2 and J3. 

Let us discuss the M = (2/3) M s magnetization plateau problem based on the above 
picture by use of the degenerate perturbation theory [T7| . The lowest energy state of 
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the jth 3-spin cluster with S^" = 1/2 is given by 

= (I hj-ihjUj+i) - 2 | hj-iUjhj+i) + I Uj-ihjhj+i)) (15) 

which is nothing but equation (@J). Here we use the symbol ipxj for convenience. The 
energy of this state is 

E 1 = -J, - l -H. (16) 

On the other hand, the state of the jth 3-spin cluster with S^} t z = 3/2 is 

ip2,j = I hj-ihjhj+i) (17) 
having the energy 

E 2 = ~J t - -H. (18) 

2 2 y 1 

Using the pseudo-spin operator Tj, where Tj 2 = Tj(Tj + l) with Tj = 1/2, we express the 

ipi,j and ip2,j states by the Tj 2 = —1/2 and = +1/2 states, respectively. We neglect 

the other 6 states of the jth 3-spin cluster, which can be justified near J 2 = J 3 = 0. 

We note that E\ = E 2 when the magnetic field is given by H = Hf^ = (3/2) Jx. The 

lowest order perturbation calculation with respect to J 2 and J 3 leads to the effective 

Hamiltonian for the pseudo-spin operator Tj 

n cS = £Ri(?i + rpj. , ) + j&fTf+j - h cS y: t> (19) 

Ji = (20) 

H«, = H - ff<°> - 5 Jl±p*L. (22) 

The ground state of the effective Hamiltonian (|19|) is either the Neel state or the SF 
state depending on whether A e s > 1 or < A c fj < 1, where A c ff is 

Acff = W\ - 6|J 2 -4J 3 |- (23) 
Thus, the Neel ground state is realized when 

-$-<*<!. (24) 
32 J 2 16 v ; 

The Neel and SF ground states in the T-picture correspond to the M = (2/3) M s 
plateauful and plateauless states in the original S-picture, respectively. The transition 
between the plateauful and plateauless states is of the BKT type |2THl2~T] . The magnetic 
field H2/3 given by the solution of iJ e fr = 0, i.e., 

„ „ (0) , 5J 2 + 22J 3 54Ji + 5J 2 + 22 J 3 

H 2/ s = H 2/3 + 36 = gg • (25) 



where 
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represents the magnetic field corresponding to the center of the magnetization plateau 
in the plateauful case, and that corresponding to the M = M s /3 magnetization in the 
plateauless case. This degenerate perturbation calculation has already been developed 
in [3[7j. The physical picture of the M = (2/3)M plateau is shown in figure ITU 



Figure 14. Physical picture of the M — (2/3)M s plateau state for the S — 1/2 DD 
chain model. Ellipses with open and closed triangles denote the 3-spin clusters with 
the ipi an d ip2 states, respectively. 

For the purpose of drawing the M = (2/3)M s plateau phase diagram, we have 
performed the numerical diagonalization of the finite size Hamiltonian by the Lanczos 
method. We can apply the LS method 012111 123 EH to this kind of plateauful- 
plateauless transition of the BKT type jHHl EH EE EH] accompanied by the SBTS. 
The plateauful-plateauless transition point can be obtained from the crossing point 
of AE (N,M = (2/3)M s ) and AE ±1 (N,M = (2/3)M s ) as functions of the quantum 
parameters. Here, AEq(N, M) and AE±i(N,M) are defined, respectively, by 

AE (N,M) = E 1 (N,M)-E (N,M) (26) 

AE ±1 (N, M) = - {E (N, M + 1) + E (N, M - 1)} - E (N, M) (27) 

where E (N,M) and Ex(N,M) are, respectively, the lowest and second lowest 
energies in the subspace of the magnetization M for the finite size system 
with iV spins. It is noted that the M = (2/3) M s state is either plateauful 
(the Neel state in the T-picture) or plateauless (the SF state in the T- 
picture) depending on whether AE (N,M = (2/3)M s ) < AE ±1 (N,M = (2/3)M s ) or 
AE (N,M = (2/3)M s ) > AE ±1 (N,M = (2/3)M s ). 

Figure [El shows the behaviors of AE (N,M = (2/3)M s ) and AE ±1 (N,M = 
(2/3)M s ) with = 18 (M s = 9) as functions of J 3 when J 2 = 0.80. The value of J 3 at 
each crossing point is expected to be an approximate value for the plateauful-plateauless 
transition point J3 in the N — > 00 limit. We have carried out, to estimate J3 , an 
extrapolation to this limit similar to that in the case of the estimation of for the 
M = M s /3 plateau phase diagram, discussed in §3, and have obtained the M = (2/3) M s 
plateau phase diagram depicted in figure HH1 In this figure a good agreement is found 
between the numerical result and the analytical prediction given by equation near 
{J21 J3) = (0, 0). We note that the point AE = ( J 3 ~ 0.25) corresponds to the Ising 
case (J^j = 0) of the pseudo-spin Hamiltonian (|T9*)l. The double degeneracy in the Neel 
ground state is realized only in the iV — > 00 limit in usual cases, whereas it is realized 
even in finite size systems for the Ising case. Thus, the point AE = 0, which means 
the doubly degenerate ground state, corresponds to the Ising case. 
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Figure 15. AE (N,M = (2/3)M s ) (closed circles) and AE ±1 (N,M = (2/3)M s ) 
(open circles) with N = 18 (M s = 9) as functions of J3 when J2 = 0.80. Approximated 
values for the plateauful-plateauless phase boundary points in the N — > 00 limit are 
known from the crossings of AE (N, M = (2/3)M s ) and AE ±1 (N, M = (2/3)M s ). The 
point where AEq(N,M — (2/3)M s ) vanishes (J3 ~ 0.25; see the inset) corresponds to 
the Ising case of the pseudo-spin Hamiltonian 119fl . 
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Figure 16. M = (2/3)Af s plateau phase diagram of the 5 = 1/2 DD chain model. 
Closed circles denote the numerical results, and the dashed lines the results of the 
degenerate perturbation calculation given by equation l|24|) . 



The expectation value (Sf) of the (2/3)M s plateau state can be obtained from the 
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above-mentioned physical picture (figure [HJ); equations (|T5j) and ()17|) yields for (S'f) 
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where three spins in a parenthesis form a 3-spin cluster. Figure El shows the behaviors 
of (Sf) which are obtained by using the DMRG method for the N = 96 system in the 
(hi J3) = (0.80, 0.27) and (J 2 , J 3 ) = (0.80, 0.30) cases. As can be readily seen, the 
behavior of (Sf) is quite consistent with equation (J2EJ). Of course, the ideal behavior 
given by this equation is expected in the A e g = 00 limit (the Ising limit) in the effective 
Hamiltonian (|19|h even if the mapping onto the T-picture is justified. In the case where 
A e ff is finite, the hybridization between the ip\ and -02 states occurs, which well explains 
the fact that (S3 -) < (S'| J _ 1 ) = (S^ +1 ) in the -0 2 -rich cluster. The behavior of (Sf) of 
the -0 2 -rich cluster in the (J 2 , J3) = (0.80, 0.27) case is nearer to the ideal behavior 
than that in the (J 2 , J 3 ) = (0.80, 0.30) case. This suggests that the parameter set in 
the former case is nearer to the Ising limit of the effective Hamiltonian (|19|) than that 
in the latter case. In fact, the Ising limit is realized when J 3 ~ 0.25 for J 2 = 0.55, as 
can be seen from figure 15. 
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Figure 17. Behavior of (Sf ) of the TV = 96 system in the M = (2/3) M a plateau state, 
where (a) (J 2 , J 3 ) = (0.80, 0.27) and (b) (J 2 , J 3 ) = (0.80, 0.30) 



5. Magnetization curves and magnetization phase diagrams 

In order to calculate the ground state magnetization curve of the present S — 1/2 
DD chain model, we have employed the DMRG method jHSEH]- The procedure of our 
DMRG calculation is briefly summarized in reference [3]. 
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We have computed E (N, M) for finite size systems with up to N = 96 (N = 216) 
spins for general (special) values of J 2 and J 3 . Once the values of E (N,M) for all 
nonnegative M's (M = 0, 1, 2, • • •, N/2) are known, the ground state magnetization 
curve can be obtained by plotting as a function of H the average magnetization 
(m) per one spin, where (m) is the value of m which yields the minimum among 
{E (N,mN)/N - #m}'s with m = 0, 1/N, 2/N, ■ ■ ■, 1/2. As is naturally expected, 
the resulting magnetization curve is a stepwisely increasing function of H. Following 
Bonner and Fisher's pioneering work [12], we may obtain, except for plateau regions, 
a satisfactorily good approximation to the ground state magnetization curve in the 
iV — > oo limit by drawing a smooth curve through the midpoints of the steps in the 
finite size results. As for the plateau regions, on the other hand, we can estimate 
the lowest and highest values of H in the N — > oo limit, denoted by H p \ and #p,h, 
respectively, giving the M = pM s (p = 0, 1/3, 2/3 or 1) plateau by extrapolating the 
corresponding finite size values to this limit. 

The magnetization curves in the N — > oo limit thus obtained for the (J 2 , J 3 ) = 
(0.80, 0.27) case and the ( J 2 , J 3 ) = (0.80, 0.55) case are shown in figure [THl It is 
noted that # p ,i and H p ^ in both cases are estimated by fitting the finite size results 
for iV = 48, 72 and 96 to a quadratic function of 1/N 2 , which is similar to equation 
(JTUJ). The H = ground state phase in the former case is the SF phase (see figure EJ), 
and therefore we have no M = plateau. Since that in the latter case is the D phase, 
on the other hand, we have in principle a finite M = plateau, although the M = 
plateau width is so narrow that we cannot see it clearly (see figure lT97 b) below). The 
wide M = M s /3 plateaux are found in both cases, while the M = (2/3) M s plateau is 
observed only in the former case, which is consistent with the M = (2/3) M s plateau 
phase diagram shown in figure El 

We have also obtained the magnetization phase diagrams on the (J3, H) plane in 
the case of J2 = 0.80 and on the ( J2, H) plane in the case of J3 = 0.55. The results are 
depicted in figure QH where H 0ih /H s , H 1/3:i /H s , H 1/3>h /H s , H 2 / 3 ,i/H s , H 2 / 3 ,h/H s as well 
as H/H s = 1.0, where H s (= Hi t \) is the saturation field, are plotted. When we estimate 
these values for general values of J 2 and J 3 , we have employed the same method as that 
for the estimations of in §3 and J 3 CT ^ in §4. It should be emphasized, however, 
that in estimating #1/3,1 and #i/3,h for special values of J2 = 1.190, 1.191, 1.192, 1.194 
and 1.195 with J 3 = 0.55, we have used the finite size results for iV = 132, 168 and 
216 instead of those for iV = 48, 72 and 96, since the iV dependence of the finite size 
results are rather strong for these J 2 and J 3 . In the case of J 2 = 0.80 the change of the 
mechanism of the M = M s /3 plateau is not observed, whereas in the case of J 3 = 0.55, it 
is clearly observed. This is quite consistent with the M = M s /3 plateau phase diagram 
shown in figure ITBI We note that in the case of J 3 = 0.55, the M = M s /3 plateau starts 
from H = when < J2 < 0.561, since the H = ground state phase is the FRI phase. 
Finally, we mention that the M = M s /3 plateau width Wnm plotted in figure fT2l is 
defined by 

W(l/3) = #i/ 3 ,h - #1/3,1- (29) 
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Figure 18. Magnetization curves in the N — > oo limit for (a) (Ja, J3) = (0.80,0.27) 
and (b) (J2, J3) = (0.80,0.55). Note that H s denotes the saturation field. The wide 
plateau at M = M s /3 is observed in both (a) and (b). The M — (2/3)M s plateau can 
be seen in (a), the parameters for which lies in the plateau region shown in figure UTtI 

6. Discussion and Summary 

We have investigated the magnetic properties of the S = 1/2 DD chain model at 
T = by use of the degenerate perturbation theory, the LS analysis of the numerical 
diagonalization data obtained by the Lanczos method, and the DMRG calculation. We 
have clarified the mechanism for the M = M s /3 and M = (2/3) M„ plateaux, and 
have also drawn the precise plateau phase diagrams by carrying out the LS analysis. 
These phase diagrams agree with the results obtained by physical considerations and 
the degenerate perturbation theory. The expectation value (Sf) of each spin estimated 
by the DMRG calculation strongly supports the plateau formation mechanisms. The 
magnetization curves for a few sets of the parameters J2 and J3 have been obtained 
by the DMRG calculation, the behavior of which are consistent with the plateau phase 
diagrams. We have also drawn the magnetization phase diagram on the (J3, H) plane 
for J 2 = 0.80 and that on the (J2, H) plane for J 3 = 0.55. 

As noted in §1, some substances are thought to be well modeled by the DD chain. 
A 3 Cu3(P04)4 with A = Ca, Sr were investigated by Drillon et al. (HIHD] and Ajiro et 
al. [TTj. The measured magnetic susceptibility x(T) [3110] showed the ferrimagnetic 
behavior. In the high field magnetization process measurement by Ajiro et al. [llj, 
the magnetization rapidly reached the wide M = M s /3 plateau as the magnetic field 
was increased. The M = M s /3 plateau was kept even at the high field limit 40T in 
their experiment. Ajiro et al. JT] performed the neutron scattering experiment for 
Sr 3 Cu3(P04)4 at T = 4K and T = 70 111K. The transition temperature of the three- 
dimensional ordering to the antiferromagnetic state of Sr 3 Cu3(P04)4 is Tn = 0.96 K. 
Although their neutron diffraction pattern at T = 70 mK was consistent with the 




Figure 19. Magnetization phase diagrams (a) on the H/H s versus J3 plane when 
J2 = 0.80 and (b) on the H/H s versus J2 plane when J3 = 0.55. In the regions A, Bx 
(X = A,B), C and D, the average magnetization (m) per one spin is equal to (the 
M = plateau region), 1/6 (the M = M s /3 plateau region), 1/3 (the M = (2/3) M a 
plateau region, and 1 (the M = M s plateau, or the saturated magnetization, region, 
respectively; the Ba and Bb regions are, respectively, for the plateaux A and B. In 
the other regions, (m) varies continuously with the increase of H. For (b), we can see 
the closing of the M = M s /3 plateau at J2 = 1.193, which corresponds to the phase 
transition between the plateau A and B phases, as expected from the M = M s /3 
phase diagram shown in figure ITU For (a), such a transition is not seen, which is also 
consistent with the M = M s /3 phase diagram. 



antiferromagnetic ordering, its magnitude was strongly reduced to almost 1/10 of that 
in usual cases. This fact suggests that the magnitude of spins on the average is about 
1/vIO ~ 1/3 of that for independent spins, which is consistent with the picture of the 
3-spin cluster formation with Sfot = ±1/2. Thus, we think, the parameter set of this 
substance lies in the FRI region on the phase diagram figure El 

Sakurai et al. [12] investigated the magnetic susceptibility, magnetization process 
(up to 28T) and 51 V NMR of Bi 4 Cu3V20i4. Unfortunately, they could not obtain 
any conclusion on the existence of the M = M s /3 plateau, because M ~ 0.27M S at 
H = 28 T. This suggests that the ground state of this substance above the three 
dimensional antiferromagnetic ordering temperature (Tn = 6 K) is the SF state. The 
ferrimagnetic behavior was not seen in the susceptibility above Tn, which is consistent 
with the SF ground state. 

Kikuchi et al. [T3Tll4|ll5j have investigated experimentally the magnetic and thermal 
properties of azurite Cu3(OH) 2 (C03)2, which seems to be modeled by the S = 1/2 DD 
chain model. Their result [13] shows that the H = ground state of azurite is the SF 
state. Furthermore, their recent result [T3] for the low temperature magnetization curve 
obtained by applying the magnetic field along the b axis (the chain axis) demonstrates 
that the M = M s /3 plateau appears in the field range of 0.5H S < H < 0.8H S with 
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H s ~ 32. 5T. These two results may imply that J 3 is considerably smaller than unity and 
J2 is considerably larger than unity, as can be seen from figure 121 since the M = M s /3 
plateau width should be much wider if J 3 ~ 1 and J 2 > 2. It is noted that the fact that 
J 3 is smaller than unity is not inconsistent with the fact [41 J that the length of the bond 
between the (3j — l)th and (3j)th sites is almost equal to that of the bond between the 
(3j)th and (3j + l)th sites, since the superexchange paths connecting S^-i and S 3 j and 
that connecting S 3 j and Sy+i are different with each other. Thus, we conjecture that 
the M = M s /3 plateau observed in azurite is the plateau B shown in figure Mjo) (see 
figure ITTj). In addition to the above experimental results, Kikuchi et al. ^2] have also 
found a plateau-like behavior at M — (2/3)M s in the dM/dH versus H curve. This 
may suggest that the value of J 2 is not so much large compared with unity (see figure 
fT6|) . Furthermore, they pB] have observed the anisotropy concerning the direction of the 
applied magnetic field in the low temperature magnetization curves. Thus, in order to 
quantitatively explain the magnetic properties of azurite, more detailed investigations, 
which take at least the antisymmetric Dzyaloshinsky-Moriya interactions [4*2714*3*] into 
account, are definitely required. We, however, believe that we have succeeded in 
explaining them at least qualitatively. 

The specific heat measurement of azurite has also been done by Kikuchi et al. |15j . 
They have found two peaks in the specific heat C{T) versus temperature T curve, which 
suggests the existence of two characteristic energies. We consider that the larger one is 
J 2 and the other is (Ji — J 3 ) 2 /2J 2 which appears in equation (j""§. 

Acknowledgements 

We would like to express our appreciation to H. Kikuchi for informing us of their 
experimental results for azurite prior to publication and for invaluable discussions. We 
are deeply grateful to Y. Ajiro for information and useful discussions on A 3 Cu3(P0 4 )4 
with A = Ca, Sr. We also thank T. Hikihara by whom the DMRG program employed 
in this study is coded, and H. Nishimori for the numerical diagonalization program 
package TITPACK Ver.2. Thanks are further due to H.-J. Mikeska, K. Takano, 
K. Kubo, H. Tanaka, H. Ohta, K. Nomura, A. Kitazawa and A. Honecker for stimulating 
discussions. This work has been partly supported by a Grant-in-Aid for Scientific 
Research on Priority Areas (B) ( "Field- Induced New Quantum Phenomena in Magnetic 
Systems") and a Grant-in-Aid for Scientific Research (C) (No. 14540329) from the 
Ministry of Education, Culture, Sports, Science and Technology of Japan. Finally, 
we thank the Supercomputer Center, Institute for Solid State Physics, University of 
Tokyo, the Information Synergy Center, Tohoku University and the Computer Room, 
Yukawa Institute for Theoretical Physics, Kyoto University for computational facilities. 

References 

[1] Ishii M, Tanaka H, Mori M, Uekusa H, Ohashi Y, Tatani K, Narumi Y and Kindo K 200 J. Phys. 
Soc. Jpn. 69 340 



Magnetic properties of the S = 1/2 distorted diamond chain atT = 



20 



Okamoto K, Tonegawa T, Takahashi Y and Kaburagi M 1999 J. Phys.:Condens. Matter 11 10485 
Takano K, Kubo K and Sakamoto H 1996 J. Phys.:Condens. Matter 8 6405 

Tonegawa T, Okamoto K, Hikihara T, Takahashi Y and Kaburagi M 2000 J. Phys. Soc. Jpn 69 
Suppl. A, 332 

Tonegawa T, Okamoto K, Hikihara T, Takahashi Y and Kaburagi M 2001 J. Phys. Chem. Solids 
62 125 

Sano K and Takano K 2000 J. Phys. Soc. Jpn 69 2710 
Honecker A and Lauchli A 2001 Phys. Rev. B 63 174407 

Fujisawa M, Yamaura J, Tanaka H, Kageyama H, Narumi Y and Kindo K 2002 J. Phys. Soc. Jpn. 
72 694 

Drillon M, Coronado E, Belaiche M and Carlin R L 1988 J. Appl. Phys. 63 3551 
Drillon M, Belaiche M, Legoll P, Aride J, Boukhari and Moqine A 1993 J. Mag. Mag. Mat. 128 
83 

Ajiro Y, Asano T, Nakaya K, Mekata M, Ohoyama K, Yamaguchi Y, Koike Y, Morii Y, Kamishima 

K, Aruga-Katori H and Goto T 2001 J. Phys. Soc. Jpn. Suppl. A, 186 
Sakurai H, Yoshimura K, Kosuge K, Tsujii N, Abe H, Kitazawa H, Kido G, Michor H and Hilscher 

G J. Phys. Soc. Jpn. 71 1161 
Kikuchi H, Fujii Y, Chiba M, Mitsudo S and Idehara T 2002 23th Int. Conf. Low Temp. Phys. 

(Aug, 2002, Hiroshima, Japan, 21CP5), to be published in Physica B 
Kamikawa T, Okubo S, Takashi K, Ohta H, Inagaki Y, Kikuchi H, Saito T, Azuma M and Takano M 

2002 23th Int. Conf. Low Temp. Phys. (Aug, 2002, Hiroshima, Japan, 21CP53), to be published 

in Physica B 

Kikuchi H and Mitsudo S, private communications 
Lieb E and Mattis D 1962 J. Math. Phys. 3 749 
Totsuka K 1998 Phys Rev. B 57 3435 

Tonegawa T and Harada I 1987 J. Phys. Soc. Jpn. 56 2153 
Okamoto K and Nomura K 1992 Phys. Lett A and references therein. 
Berezinskii Z L 1971 Zh. Eksp. Teor. Fiz. 61 1144 (Sov. Phys.-JETP 34 610). 
Kosterlitz J M and Thouless D J 1973 J. Phys. C: Solid State Phys. 6 1181 
Okamoto K and Kitazawa A 1999 J. Phys. A: Math. Gen. 32 4601 
Honecker A 1999 Phys. Rev. B 59 6790 

Oshikawa M, Yamanaka M and Affleck I 1997 Phys. Rev. Lett. 78 1984 
Okamoto K 1992 Solid State Commun. 83 1098 
Okamoto K 2002 Prog. Theor. Phys. Suppl. No. 145, 113 

Nomura K and Kitazawa A 2002 Proceedings of French- Japanese Symposium on Quantum 
Properties of Low- Dimensional Antiferromagnets (eds. Y. Ajiro and J. -P. Boucher, Kyushu 
University press, 2002): |cond-m at/02010T2| 

Kitazawa A 1997 J. Phys. A: math. Gen. 30 L285 

Kitazawa A and Okamoto K 1999 J. Phys.: Cond. Matter 11 9765 

Kitazawa A and Okamoto K 2000 Phys. Rev. B 62 940 

Sakai T and Okamoto K 2002 Phys. Rev. B 65 214403 

White S R 1992 Phys. Rev. Lett. 69 2863 

White S R 1993 Phys. Rev. B 48 10345 

Nomura K and Okamoto K 1994 J. Phys. A: Math. Gen. 27 5773 
Nomura K 1995 J. Phys. A: Math. Gen. 2 5451 

Okazaki N, Okamoto K and Sakai T 2000 J. Phys. Soc. Jpn. 69 2419 
Okazaki N, Okamoto K and Sakai T 2001 J. Phys. Soc. Jpn. 70 636 

Nakasu A, Totsuka K, Hasegawa Y, Okamoto K and Sakai T 2001 J. Phys.: Condens. Matter 13 
7421 

Okamoto K, Okazaki N and Sakai T 2002 J. Phys. Soc. Jpn. 72 Suppl. 196 
Bonner J C and Fisher M E 1964 Phys. Rev. 135 A640 



Magnetic properties of the S = 1/2 distorted diamond chain at T 



[41] Zigan F and Schuster H D 1972 Z. Krist. 135 416 
[42] Dzyaloshinsky I 1958 J. Phys. Chem. Solids 4 241 
[43] Moriya T 1960 Phys. Rev. 120 91 



